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ABSTRACT 

In  the  first  part  of  this  paper  we  study  in  detail  the  properties  of  the 
divergence  operator  acting  on  continuous  piecewise  polynomials;  more 
specifically,  we  characterize  the  range  and  prove  the  existence  of  a  maximal 
right-inverse  whose  norm  grows  at  most  algebraically  with  the  degree  of  the 
piecewise  polynomials. 

In  the  last  part  of  this  paper  we  apply  these  results  to  the  p-version  of 
the  Finite  Element  Method  for  a  nearly  incompressible  material  with 
homogeneous  Dirichlet  boundary  conditions.  We  show  that  the  p-version 
maintains  optimal  convergence  rates  in  the  limit  as  the  Poisson  ratio 
approaches  V2  •  This  fact  eliminates  the  need  for  any  "reduced  integration" 
such  as  customarily  used  in  connection  with  the  more  standard  h-version  of  the 
Finite  Element  Method. 
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A  RIGHT-INVERSE  FOR  THE  DIVERGENCE  OPERATOR  IN 
SPACES  OF  PIECEWISE  POLYNOMIALS.  APPLICATION 
TO  THE  p- VERSION  OF  THE  FINITE  ELEMENT  METHOD. 

Michael  Vogelius* 

1.  Introduction 

For  many  finite  element  discretizations  of  p.D.E. 's 
estimates  of  the  error  rely  on  a  detailed  study  of  the  divergence 
operator  on  corresponding  sets  of  piecewise  polynomials ‘Examples 
of  this  are  mixed  formulations  for  the  Laplace  equation  or  the 
equations  of  elasticity  (cf.  [6],  19])  and  discretizations  of  the 
Stokes  problem  (cf .  [10] )  . 

Finite  element  methods  in  general  are  based  on  some  partition 
of  the  domain  and  convergent  solutions  are  typically  obtained  by 
letting  the  mesh  size  h  tend  to  0  using  piecewise  polynomials  of 
some  fixed  degree  (the  h-version) .  Recently  a  different  version 
has  received  quite  a  bit  of  attention  (cf.  11,3,4,7]),  namely  the  so-called 
p-version  of  the  Finite  Element  Method.  The  p-version  is  based  on  a 
fixed  triangulation  and  the  degree  of  the  piecewise  polynomials  is 
increased  in  order  to  give  convergence. 

In  section  2  of  this  paper  we  analyze  the  divergence  operator 
with  special  emphasis  on  applications  for  the  p-version,  to  be 
more  precise 

(i)  we  characterize  the  range  of  the  divergence  operator 
acting  on  fields  that  are  continuous  piecewise  poly¬ 
nomials  of  degree  £  p+1  on  some  triangulation  of  a 
domain  fi  (for  p  ^  3)  . 
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(ii)  we  show  that  a  maximal  right-inverse  exists,  which 
takes  values  in  continuous  piecewise  polynomial 
fields  of  degree  _<  p+1  ,  and  which  has  an  operator 
norm  that  grows  at  most  algebraically  with  p. 

The  first  of  these  two  results  is  by  far  the  easiest,  it 
may  be  obtained  from  a  purely  combinatorial  argument  and  use  of 
the  Grassmann  dimension  formula.  Since  (i)  follows  very  naturally 
as  an  intermediate  result  in  our  proof  of  the  second  (and  much 
harder)  assertion  we  shall  only  give  a  sketch  of  the  selfcontained 
combinatorial  argument  (cf.  Remark  2.1). 

The  reference  [9]  contains  an  explicit  construction  of  a  right- 
inverse  for  the  divergence  operator  for  piecewise  polynomials  of 
arbitrarily  high  degree,  but  the  fields  defined  that  way  are  not 
continuous  across  inter-element  boundaries,  only  the  normal 
components  are. 

In  section  3  we  apply  the  results  (i)  and  (ii)  to  a  particular 
discretization  by  the  p-version.  The  continuous  model  is  that  of 
a  nearly  incompressible  linearly  elastic  material  (Poisson  ratio 
v  <v.  i )  with  boundary  conditions  that  are  homogeneous  Dirichlet. 

We  show  that  the  discrete  solution  converges  at  an  optimal  rate 
even  when  v  is  extremely  close  to  j  »  this  is  in  sharp  contrast  to 
the  h-version  where  one  encounters  significantly  reduced  convergence 
rates  for  nearly  incompressible  materials  (we  refer  to  [11]  for  a 
more  detailed  discussion) .  In  [11]  we  relied  on  the  results  (i) 
and  (ii)  to  analyze  the  similar  problem  for  the  case  of  natural 
boundary  conditions  on  a  smooth  domain. 

Let  D  be  a  bounded  polygonal  domain  in  Rn,  n  =  1,2.  . 

Throughout  this  paper  H  (D)  denotes  the  standard  Sobolev  space 

2 

of  functions  that  have  all  derivatives  of  order  <  k  in  I.  (D) . 


Whenever  no  ambiguity  is  possible  the  same  notation  shall  also 
be  used  for  vector  valued  functions  whose  components  are  in 
H  (D) .  Spaces  with  noninteger  exponents  are  defined  by  complex 
interpolation  {cf.  (5]).  The  norm  on  the  space  HS(D)  is 

O 

denoted  1*1  „  and  HS(D)  refers  to  the  closure  of  cT(D) 

s »  d  u 

in  HS(D) . 
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2 .  A  maximal  right-inverse  for  the  divergence  operator. 

2 

Let  n  be  a  bounded  polygonal  domain  in  JR  and  let  £  denote 
a  triangulation  of  £2  .  We  assume  that  £  has  the  property  that  no 
vertex  of  a  triangle  lies  in  the  interior  of  a  side  of  another. 

Before  treating  piecewise  polynomials  on  the  triangulation  £ 
we  establish  a  few  preliminary  facts  concerning  polynomials  on  a 
single  triangle.  In  what  follows  T0  shall  always  denote  the 
reference  triangle  ((x,y)|o  <  x,  0  <  y,  x+y  <  1}  .  |jj-  is  the 
derivative  along  the  outward  normal  to  3Tg  and  is  the  derivative 
tangential  to  and  counterclockwise. 


Lemma  2 ■ 1 

Let  Qp+*(x)  and  rP (x)  be  two  polynomials  in  x  of  degree  <  p+1 
and  <_  p  respectively.  Assume  that  0  and  1  are  zeros  of  multiplicity 
at  least  3  for  Qp+1  and  at  least  2  for  RP. 

There  exists  a  polynomial  Sp+1(x,y)  in  x  and  y  of  degree 


=  0  on  the  two  sides 

(x)  , 

'P  f  v  1  nn  t- Hn  ei  v±h 


<  p+1 

satisfying 

(i) 

sptl  -  Ih  s"*1  ■ 

x  =  0  ,  y  =  0, 

(ii) 

Sp+1(x,y)  =  Qp+1 

In  sP+1<*'y>  = 

(iii)  «sp+1l  <c<p+uK  (|QP+I'0f[0fll  * 


,rP,<M0>1]) 

with  constants  C  and  K  that  are  independent  of  p, 
qP4-1  and  RP  . 
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Proof:  Since  0  and  1  are  zeros  of  multiplicity  j>  2  for  both 

qp+*  and  Rp  we  can  write 

QP+1(x)  =  x2(l-x)2  QP_3(x) 

Rp{x)  =  x2(l-x)2  Rp"4(x)  . 


Define  Tp+^  by 

Tp+1(x,y)  =  x2y 2  QP_3(x)  . 

Then 

(1)  TP+1  _  |_  tp+1  =  0  on  the 

<?n 

two  sides  x  =  0  ,  y  =  0  . 

(2)  Tp+1 (x,y)  =  Qp+1 (x)  on 
the  side  x+y  =  1  . 

By  integration  we  get 


jf  (Tp+1(x,y) ) 2  dx  dy  =  j  x4(Qp_3(x))2  j  y4  dy  dx 
'  T0  0  0 
1  1 
=  |  j  x4(l-x)5(QP'3(x))2dx  <  |  |  (QP+1 (x) ) 2  dx  , 


x .  e. 


<3) 


,tP+1|o,  r0  i  c 


Due  to  the  fact  that  0  and  1  are  zeros  of  multiplicity  at 
least  3  for  Qp+*,  we  can  also  write 
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h  TP+1(X>y)  y=l-x  =  x2(l-^)2  ^(x)  . 


Define  Up+^  by 


UP+1(x,y)  =  /I  x2y2  (x+y-1) (Rp_4 (x)  -  Tp-4 (x) ) 


Then 

(4)  Up+^  =  Up+^  =0  on  the 
two  sides  x  =  0,  y  =  0  . 

(5)  Up+1  =  o  and 

Up+1(x,y)  =  Rp(x)  -  Tp+1(x,y) 

on  the  side  x+y  =  1. 


By  integration  as  before, 


,P+1. 


o,  t0  <  c(,rP,o,(o,ii  +  >Ik  tP+1L=1_x  * 0 ,  [0 , 1  ] 5 


The  standard  trace  theorem  therefore  gives 

(6)  IUP+^I  <  C  ( | rp |  +  .mP+l,  ) 

°'T0  -  0,[0,1]  +  •*  't.T  ' 


for  any  t  >  3/2  . 

It  is  known,  and  used  repeatedly  in  this  paper,  that 

ITP+1I.  T  <  C.  (p+l)2t  I Tp+1 I „  j 
-  t  0,Tq 

■  for  any  polynomial  Tp+*  of  degree  <  p+1,  and  any  t  >  0. 

(The  reference  [4]  contains  a  proof  of  this  inequality  on  a 


square.  By  a  localization  argument  one  can  extend  to  other 
type  domains,  e.g.  triangular).  Inserting  into  (6)  and  using  (3) 
we  thus  obtain 


V  T0  -  ctlp+1|2t('RPY [o,i]  *  '^Vnu' 


for  any  t  >  3/2. 

Finally  let 

Sp+1  =  Tp+1  +  U^1  . 

It  is  then  clear  from  (1),(2),(4)  and  (5)  that  Sp+1  satisfies 
the  requirements  (i)  and  (ii)  .  The  estimate  (iii)  follows  from 
(3)  and  (7) .As  K  we  can  use  any  number  >3.  □ 


Lemma  2.2.  Let  Cp  (x)  and  Rp(x)  be  two  polynomials  in  x  of  degree 
—  p.  Assume  that  0  and  1  are  zeros  of  multiplicity  at  least  2  for 
Qp  and  Rp.  Assume  furthermore  that  [  Qp(x)  dx  =  0. 

Jo 

There  exists  a  polynomial  Sp+1(x,v)  in  x  and  y  of  degree 
<_  p+1,  satisfying 


■sP+  'o  T  £  C(p+1)K  (»Qpl 
'  0 


0, [0,1] 


+  IRPI 


0, [0 , 1 J 


with  constants  C  and  K  that  are  independent 
of  p ,  Qp  and  Rp  . 
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Proof:  Define  Qp+1(x)  =  -  [  Qp(s)ds  and  construct  Sp+1tx,y) 
according  to  the  previous  lemma  with  boundary  data  Qp  and  Rp. 
This  Sp+1  has  the  desired  properties.  □ 


lemma  2.3.  Let  T  be  a  triangle  two  sides  of  which  lie  on  the  lines 
a^x  +  b^y  +c^=0,  i  =  1,2.  Let  cP  be  a  polynomial  of  the  form 


Qp(x,y)  =  (axx  +  bjy  +  Cj^)  (a2x  +  b2y  +  c2)  Rp  (x,y)  , 

where  Rp-^  is  a  polynomial  in  x  and  y  of  degree  <  p-3.  Assume 
furthermore  that 


JJ  QP(x,y)dx  dy 

There  exists  a  field 

degree 

<  p+1# satisfying 

(i) 

Up+1  =0  on  JT 

(ii) 

V-I2P+1  =  QP  in  T 

(iii) 

iup+1i1>t  <  C (p+1 > K 

with  constants  C  and  K  that  are  independent  of 
p  and  Qp. 


9 


Proof :  By  an  affine  transformation  T  can  be  mapped  onto  Tg  in 
such  a  way  that  a^  +  b1y  +  c^  =  0  goes  to  x  =  0  and  a2x  +  b2y 
+  c2  =  0  to  y  =  0.  It  therefore  suffices  to  prove  the  result  on 
Tg  with  Qp  given  by 

QP(x,y)  =  xy2  Rp_3(x,y)  . 


Define  w 


P+1 


as  follows 


Wp+3 (x ,y)  =  |  Qp(s,y)ds 

0 

,P+1, 


W' 


(x,y)  =  0  . 


This  ensures  that 

HP+! 


(x,y)  =  x2y2  Wp_3(x,y) 


•wp+1  =  op 


(8) 


v-w1 

wP+1 


(T  and 


1 KP+ 1 • ,  T  £  C(p+1)2  !QpH 

1  •  1  n  U  r  I  r 


Wp+1  is  therefore  0  on  the  two  sides  x  =  0  and  y  =  0.  To  make  it 
vanish  also  along  the  side  x+y  =  1  we  now  define  an  appropriate 
correction  term.  Let  Sp+2  denote  the  polynomial  constructed  in 
Lemma  2.2  with  boundary  data  n-wp+1 (x ,  1-x) ,  -s-Wp+1 (x ,1-x)  for  the 
tangential  derivative  and  the  normal  derivative  respectively  (here 


we  use  the  fact  that 


Q‘  (x,y)dx  dy  =  0  )  .  Then 


UP+1(x,y)  =  Wp+1(x,y)  -  (|~  Sp+2 


■  -  k  sp*2> 
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satisfies  the  requirements  (i)  and  (ii)  ,  By  integration  it  also 
immediately  follows  that 


IW 


P+1 1 


jy=l-x'o, [0 


,1]  -C,qP'o, 


Combining  this  with  the  estimate  (iii)  of  Lemma  2.2  one  finds 


(9) 


"w  $P*2'  -  fc  sP*2  "l.T  <  C‘P*1)4  lsP*2|0,  T 


<  C(p+1)K  IQPln  T 
u,  i  o 

for  any  K  >  7.  The  estimates  (8)  and  (9)  easily  yield  (iii) . 

□ 

The  following  lemma  characterizes  the  range  of  the  divergence 
operator  for  the  case  of  a  single  triangle  and  polynomial  fields  with 
homogeneous  boundary  conditions.  In  addition  it  establishes  a  bound  on 
the  norm  of  a  maximal  right- inverse  for  V-.  As  will  be  clear 
from  the  proof, the  latter  is  the  most  difficult. 


Lemma  2.4.  Let  T  be  a  triangle  the  vertices  of  which  have 
coordinates  (x^.y^),  1  i  <3.  Let  Qp(x,v)  be  a  polynomial  in 
x  and  y  of  degree  £  p  such  that  Qp(x^,y^)  =  0,  1  £  i  £  3,  and 
||  Qp(x,y)dx  dy  =  0. 

There  exists  a  field  U_p+*  =  (Up+^,  U^+1)  of  polynomials 
of  degree  £  p+1 , satisfying 


(i) 

Up+1  =  0 

on  3  T 

(ii) 

7 -UP+1  = 

Qp  in  T 

(iii) 

1  f 

T  £C(P+DK  »Qp»0it 

with  constants  C  and  K  that 
are  independent  of  p  and  Qp. 
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Proof :  The  lemma  is  trivial  for  p  <_  1,  so  in  the  following  we 
shall  always  assume  that  p  ^  2.  Let  Pp+^  denote  the  polynomials 
in  x  and  y  of  degree  <_  p+1  and  let  Pp+^  denote  pp+^  n  ft1  (7 )  .  it 
is  easy  to  check  that  the  divergence  operator  maps  Pp+^xPp+^  into 


the  subspace  of  P* 


characterized  by  Qp (x^,y^)  =  0,  1  <_  i  <_  3,  and 


|j  Qp(x,y)dx  dy  =  0.  The  dimension  of  Pp+^x  p p+^  is  (p-l)p.  The 


null  space  of  7* 


is  for  p  ^  4  given  by 
3 


{(fy  QP+2,  -  QP+2)  |QP+2(x,y)  =  |  j  (aix+biv+ci)2  RH(x,y)  ,  RH£fp’4}, 

where  aix  +  b±y  +  c±  =  0,  1  <  i  <  3,  are  the 

equations  for  the  3  sides  of  T  .  For  p  =  2,3  the  nullspace  consists 
of  0  only.  Hence  in  general  the  null  space  has  dimension  (p— 3) (p-2)/2. 

By  using  the  Grassmann  dimension  formula  we  now  get  that  the  range  of 

V*  must  have  dimension  (p-l)p  -  (p-3)(p-2)/2  =  (p+l)(p+2)/2  -  4. 

Since  the  subspace  of  Pp  characterized  by  Qp(x^,y^)  =  0,  1  <  i  <  3, 
and  ||  Qp(x,y)dx  dy  =  0  has  dimension  exactly  (p+l)(p+2)/2  -  4 ,  we 

soace.  In  other  words, 
and 

XX  —  — 

||T  QP(x,y)dx  dy  =  0  one  can  find  Up+1,  a  field  of  polynomials  of 
degree  <_  p+1,  with  the  properties  (i)  and  (ii)  . 

The  estimate  (iii)  does  not  follow  from  the  previous  argument. 

The  existence  of  a  bounded  maximal  right-inverse  for  any  fixed  p  is 
guaranteed,  what  we  do  not  know  is  how  the  norm  depends  on  p. 

In  order  to  prove  (iii)  it  is  certainly  sufficient  to  consider 

p>5.  Let  Qp  be  a  polynomial  of  degree <_ p  which  satisfies  Qp(xi,yi)  =  0, 

1  i  i  i  3<  and  ff  (^(XfyMx  dy  =  0.  One  can  always  find  a  polynomial 


conclude  that  V-  maps  pp+1  x  pp+1  onto  this 
given  Qp  e  Pp  which  satisfies  QP(x^,y^)  =  0,  1  <  i  <  3, 


ililriTlBTm’  ■  tit  t+at*  L  ii  i— 


i 
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R5  of  degree  £  5  such  that 

(10)  <|j>j  (fc)k  (qP  -  *5) (xi»yi)  =  0 

for  0  £  j+k  £  2,  1  £  i  £  3, 

(11)  Jj  (QP-R5) (x,y)dx  dy  =  0  and 
‘  JT 

(12)  ,r5,o,t  -  c,Qp,t,r  for  any  t  >  3. 

Due  to  (10)  and  (11),  R^(x^,y^)  =0,  1  £  i  £  3,  and  jj  R^(x,y)dx  dy  =  0. 

The  combinatorial  argument  at  the  beginning  of  this  proof  therefore 
gives  the  existence  of  a  field  V6  with  the  properties 

V6  =  0  on  3T  t 

V*V6  =  R5  in  T  and 

(13)  .v'.1T  <  CIR5.0  r  . 

Based  on  (12)  and  (13)  we  get 

,y‘,1T  <clo»ltfI 

£  C(p+l)2t  IQPI0  T 

for  any  t  >  3. 

Thus  it  suffices  to  prove  the  assertion  (iii)  for  polynomials 
Qp  that  additionally  vanish  up  to  and  including  second  derivatives  at 
the  vertices  of  T.  Such  a  Qp  may,  according  to  Lemma  2.1,  be  written  as 
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Qp(x,y)  =  (aj^x  +  bxy  +  2  Ca2x  +  b2y  +  c2)  2  Q^”4(x,y) 

+  (a2x +  b2y + c2) 2 (a3x + b3y + c3) 2  Op  4(x,y) 

+  (a3x +  b3y +  c3) 2 (a3x +  b3y + c^) 2  Qp-4(x,y)  , 

2 

where  the  integral  of  each  term  over  T  vanishes,  and  the  L  -norm  of 
each  term  is  bounded  by 

C (p+ 1) K  I QP I 0  T  . 

An  application  of  Lemma  2.3  term  by  term  now  leads  to  the  desired 
result.  □ 

In  order  to  extend  the  previous  result  from  a  single  triangle 
to  the  entire  triangulation  we  shall  need  one  more  lemma.  Let  us 
first  introduce  some  notation.  p  I P 1 » 0  denotes  the  set  of  continuous 
functions  on  ft,  whose  restriction  to  each  triangle  is  given  by  a 
polynomial  of  degree  at  most  p.  By  an  internal  vertex  of  E  we 
understand  a  vertex  which  lies  in  ft  (not  on  3ft)  .  p^Pl'  ^ 
denotes  the  set  of  functions,  whose  restriction  to  each  triangle 
is  given  by  a  polynomial  of  degree  at  most  p  (no  inter-element 
continuity  requirements) ,  with  the  additional  property  that: 


R.l.  "At  any  internal  vertex  (xQ,y0)  of  the  trianqulation ,  where 
four  triangles  come  together  as  in  Fig.  1 

(i.e.,  the  sides  meeting  at  (x^y^)  fall  on  two  straight  lines  ) 
the  corresponding  values  (Qp  (xo'Yo^i  =  l  satisfv 

I  Q^(Xo,y0)  =  I  ^i*x0'v0^  ’ 

i  even  i  odd 


2 


Fig.  1 

Lemma  2.5.  Let  Vp  be  an  element  of  PIp1'-1.  There  exists 
W3  €  pf3! /0 x  p[3] ,0  such  that 

(i)  Vp-V-W3  =  0  at  all  vertices 
of  the  triangulation  L  . 

(ii)  iw3.1>fi  <  C(P+1)K  ivp.0(fi  , 

with  constants  C  and  K  that  are 
independent  of  p,  Vp  . 

Proof:  Let  us  consider  two  triangles  placed  adjacent  to  each  other 
(as  in  Fig.  2) .  Q 
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3 

V.'e  easily  construct  a  field  U  ,  given  by  third  order  poly¬ 
nomials  in  each  triangle,  continuous  in  u  T^,  vanishing  along 
3(T3  u  T21  (the  outer  boundary)  and  such  that  each  component  has 
derivative  0  in  the  direction  n2  at  the  vertex  O'  and  a  prescribed 
derivative  in  the  direction  n2  at  the  vertex  0.  If 


[? 


=  au  nx  +  a12  n2 


a21  nl  +  a22  n2  ' 


then 

(14) 


’•a3  -  -12  3nJ  “l  *  ‘22  Isj  »2 


in  T ^  at  the  vertex  0.  Similarly  if 


r° 

Lb 


-  bn  n2  +  b12  n3 


=  b2i  n2  +  b22  n3  , 


then 

(15) 


V*^3  =  bll  InT  U1  +  b21  TnZ  U2 
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*22 


and 


521 


are 


in  T2  at  the  vertex  0, 

It  is  easy  to  see  that  the  vectors 
linearly  independent  except  in  case 
n3  =  - ■  Since  we  have  the  freedom  to  prescribe  ypp-  at  the 

vertex  0,  (14)  and  (15)  enable  us  to  prescribe  any  pair  of  values 
for  V-U  (0),  one  value  in  and  another  in  T 2  ,  provided  and 
are  linearly  independent.  If  ^  and  n3  are  linearly  dependent  then 
the  formulas  (14)  and  (15)  prescribe  the  same  value  for  V>ubo) 
in  T3  as  well  as  7"2. 


Let  us  now  apply  this  preliminary  result  to  an  internal 


Fig.  3 


First  suppose  that  there  exists  a  pair  of  adjacent  triangles  such 
that  the  combined  angle  is  not  n(e.g.  tn  and  T,  in  Fig. 3  ). 

By  the  procedure  just  described,  we  may  prescribe  the  same  value 
for  the  divergence  at  0  in  the  two  triangles  and  Now  we 

turn  counterclockwise  and  do  the  same  in  T2  and  T 3.  Continuing 
that  way  until  we  reach  the  last  triangle  TN  and  performing  a 
summation,  we  obtain  the  following  set  of  prescribed  values: 

T1  !  al 

T2  !  al  +  a2 

T3  :  a2  +  a3 

TN-1  :  aN-2  +  aN-l 
TN  :  aN-l 

Since  the  a^,...,aN_3  may  be  chosen  arbitrarily,  this 
clearly  gives  the  possibility  of  assigning  any  set  of  values  in 
the  triangles  through  TN_^.  By  assumption  TN  and  do  not 
have  a  combined  angle  of  *  so  we  can  also  prescribe  a  value  a„ 
in  and  0  in  .This  permits  us  to  construct  any  set 
of  values  for  the  divergence  at  0  in  the  triangles  through  TN. 


17 


The  only  way  there  can  be  no  two  adjacent  triangles  with 
combined  angle  ^  n  is  if  the  vertex  is  the  intersection  of  two  straight 
lines  (Fig.  1).  Taking  two  triangles  at  a  time  as  before,  we  can  now 
prescribe 


i.e.,  exactly  those  sets  of  values  that  satisfy 

2  value(T^)  =  £  value(T^) 

i  even  i  odd 

This  is  precisely  the  constraint  imposed  on  »_1  at  such  vertices. 

For  any  Vp  G  p[P]'-1  we  are  therefore  able  to  find  a  field  U3 
(continuous  in  ft)  such  that  Vp  -  V*U3  =  0  at  a  particular  internal 
vertex  and  V-u3  =  0  at  all  other  vertices.  Our  construction 
has  the  property  that 

■u3iifft  <  C(P+1)K  «vp.0(n 

for  any  K  >  2. 

By  a  similar  construction  we  can  also  obtain  a  u3  such  that 
VP- V • u 3  =  0  at  a  particular  vertex  on  3ft  and  V*U3  =  0  at  all  other 
vertices.  Aiding  over  all  vertices  of  t  we  qet  the  desired  field. 

□ 

We  are  now  finally  in  a  position  to  prove  the  main  result  of 
this  section. 

Theorem  2.1.  Let  Pp  denote  the  space  V*  (  ptp+^  '  plp+l),0j. 

There  exists  a  linear  operator  L  from  Pp  into  pIP+D'0  „p[p+l]«0 


such  that 
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<i)  7.<tpVP)=VP  VvPepP 

(ii)  <tp  Vpl1>n  <  C(P+1)K  lvpl0n  V  VP  €  PP 

with  constants  C  and  K  that 
are  independent  of  p,  Vp. 

Furthermore,  there  exists  p^  such  that 

pP  =  p[pJ /-I 

for  any  p  >  pQ. 

Proof :  We  begin  by  showing  that  given  Vp  e  ptPl*  there  exists 
uP+1  e  plP+11,0  x  ptp+l],0  such  that 

(16)  lvp- v-up+1i0/n  <  c*(p+l)_t  ,  t*  >  0. 


It  is  sufficient  to  prove  (16)  for  p  >  9.  Because  of 
Lemma  2.5  it  obviously  suffices  to  consider  Vp  that  satisfy 
Vp  ■=  0  at  all  vertices  of  E  and  verify  that 

(17)  ,vP~  v*yp+1,o,n  -  c(p+i)_t  'vPIq  jj 

for  some  t  >  2.  Construct  a  piecewise  ninth  degree  polynomial 

~g  , 

V  that  vanishes  up  to  and  including  second  derivatives  on 

all  sides  of  E  (and  therefore  at  all  vertices)  with 


I/, 


(Vp  -  v9)  dx  dy  =  0 


on  every  triangle  T.This  satisfies  l\r*3  <_  ClVpl0 

Application  of  Lemma  2.4  separately  to  each  triangle 
shows  that  we  can  find  Wp+*  e  plp+lJrO  x  p[p+l],0  with 

V-.wp+I  =  Vp  -  V9 

(18)  rvM  X  „ 

•wp  1,lfn  I  c(P+i)K  ivpi0>n  . 
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It  is  possible  to  determine  0  such  that 


and 


*  v®  in  ft 

'•'Mi  ^Vo 


Hence  there  exist  polynomials  *p+2  of  degree  <_  p+2  with 

IA(4>- *p+2)»q <  CCp+l)“3  lV9l3fn 
<  c(P+i)-3  ivPio  n  . 

Therefore  i*>+1  »  V$p+2  satisfies 

(19)  ’y9-  ^.yP+1l0  R  <  c(p+l)-3  IVplon  . 

Let  Up'v  (  h<  given  by  up+1  «  Wp+1+Yp+1,  a  combination  of  (18) 
and  vl9)  then  yields  (17)  (with  t  =  3) . 

It  is  easily  checked  that  the  inclusion  Pp  C  ptPl*”1  holds. 
On  the  other  hand  let  pQ  be  such  that  C  (p0+l)  r  <  1 
(C*,  t*  are  here  the  specific  constants  in  (16))  it  follows 
then  directly  from  (16)  that  the  orthogonal  complement  of  Pp 
in  pfP^'-1  contains  0  only,  i.e.  Pp  =  P^'-1,for  p  >  pQ. 

We  now  go  back  to  (18) }  instead  of  approximating 
V9  e  V  q°  ,  qQ  -  max  {pQ,9},  we  may  at  this  point  use 

*  q0+1  iq0+1]'°  iv11'0 

the  fact  that  we  can  find  Z  e  p  xpu  such  that 

<Jn+l  -o 


and 


V  •  z 
l^*1 


*i,n  -  c,^9,o,n 


CIVPI 


o,n  • 


Defining  Up+^  =  Wp+^  +  z**®  e  ptP+H»0  x  p(p+l),0  , 


(P*.^)  we  get 
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(20)  V  •  UP+1  =  VP  , 

and 

(21)  ,yp+1,iJn  i  C(p+1)K  ivpi0>n  . 

In  the  construction  leading  to  (20)  and  (21)  we  assume  that 
VP  €  plpl*“l  satisfies  Vp  =  0  at  all  vertices  of  I.  Applying 
Lemma  2.5  we  obtain  that  (20)  and  (21)  may  be  satisfied  for 
any  Vp  e  p^PJ'“*  (p  >_  q^) .  Let  L  Vp  denote  the  field  in 
p[p+l]»0  x  ptp+1] .0  With  minimai  H1-norm  for  which  (20)  holds. 
Then  L  is  linear,  (i)  holds  by  definition  and  (ii)  follows 
from  (21) .  This  establishes  the  existence  of  an  apropriate 
operator  Lp  provided  p  >  q0,and  that  is  obviously  sufficient.  D 


Remark  2.1.  If  the  domain  is  simply  connected  then  the  divergence 
operator  acting  on  Plp+1]  '°  *  plp+11 '°  has  as  its  null  space 


^  S?  2  '  ~  fx  sP+2M  SP+2  is  a  C1  piecewise  polynomial  of  degree 
<  p+2}.  From  the  formula  for  the  dimension  of  the  space  of  C1 
piecewise  polynomials  of  degree  <  q  (q>5)  that  is  given  in  [  8j 


one  can  compute  the  dimension  of  this  null  space  in  case  p  >  3 
(the  vertices  that  fall  at  the  intersectionof  two  straight  lines  plav 
a  special  role  in  this  formula)  .  By  the  Grassmann  dimension  formula 
one  finds  that  Pp  =  fQr  any  p  >  3.  This  argument  also 

extends  to  cover  polygonal  domains  ft  that  are  not  simply 
connected, still  showing  that  Pp  =  pfPJ»“l  fDr  any  p  >  3. 
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Remark  2.2.  Let  (2  be  an  arbitrary  bounded  Lipschitz  domain. 
We  now  permit  a  triangulation  of  ft  to  have  triangles  with 
one  curved  side  (in  common  with  3(5) . 

Assume  that  any  of  the  triangles  with  one  curved  side  can  be 
mapped  affinely  onto  a  Lipschitz  domain  of  the  form 

{ (x,y)  |  0  <  x  <  1,  0  <  y  <  g(x)>  , 

where  g  is  monotone  and  satisfies  g{0)  =  1,  g(l)  =  0. 

With  this  assumption  one  may  easily  check  that  Theorem  2.1 
remains  valid  also  for  triangulations  of  nonpolygonal 
domains  (5. 

Remark  2.3.  The  proof  of  Theorem  2.1  immediately  leads  to 
a  similar  theorem  for  fields  that  satisfy  homogeneous  Dirichlet 
boundary  conditions  on  a  polygonal  domain.  The  space 
corresponding  to  *  will  now  be  defined  through  the 

two  additional  requirements: 

R.2.  "At  any  vertex  (xQ,y0)  on  3(2  where  k  triangles, 

1  <  k  <  4,  come  together  as  in  one  of  the  cases 
shown  in  Fig.  4  (i.e.  the  sides  meeting  at  (x0,yQ) 
fall  on  two  straight  lines) , 


i  <  k  i  <  k 


P  P 

l  °i(xo'yo>  =  E  Qi(xo 

i  even  1  u  0  i  odd  10 


R.3. 


JJ  QP(x,y) 
(2 


dx  dy 


0  . 


# 
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AoDlication  to  the  p-version  F.E.M. 


In  this  section  we  study  the  behavior  of  the  p-version  of  the 
Finite  Element  Method  for  the  case  of  a  nearly  incompressible 
material.  The  continuous  problem  is  that  of  plane  strain  linear 
elasticity  on  a  bounded, simply  connected  polygonal  domain  ft. 

We  furthermore  assume  that  all  corners  of  ft  have  interior  angles 
<  2it.  Formulated  in  terms  of  displacements  uv  =  (u^  uv  2) 

the  continuous  problem  is 


(22)  -Auv  -  m-uv)  *  F  in  ft, 

where  F  represents  an  external  force  (to  be  specific  f  is  the 
external  force  scaled  by  2(l+v)/E).  E  >  0  denotes  the  Young  modulus 
and  0  <  v  <  the  Poisson  ratio;  the  case  that  v  is  very  close  to  H 
corresponds  to  a  nearly  incompressible  material.  The  boundary 
conditions  are 

(2  3)  uv  =  0  on  3ft  . 

The  analysis  presented  here  could  easily  be  extended  to  the  case 
that  homogeneous  Dirichlet  data  is  given  on  part  of  the  boundary 
and  natural  (stress)  boundary  conditions  on  the  rest. 

Solving  the  boundary  value  problem  (22)  and  (23)  is  equivalent  to 
minimising  the  energy  expression 


<241  (J  [if*  V2  *  <1?  V2  *  I'll  v2  *  h  Vl’2  *  Av,V'^2]  *■ dv 


-  j|  F  v  dx  dy 
fi 

in  (Q)  .  (We  minimally  assume  that  F  e  H  the  dual  of  H  (fJ).) 

For  values  of  very  close  to  \  (22)  and  (23)  may  be  viewed  as 
a  penalized  version  of  the  Stokes  problem 

(25)'  -  AU  +  VP  =  F  and 

V*U  =0  in  n  ,  with 

U  =  0  on  • 

(It  is  well  known  that  (25)  has  a  unique  solution  U  6  H  (fi) , 

P  6  L 2(ft)  /  {constants}  provided  F  €  H  ^P),  cf.  [10].) 


We  shall  make  use  of  the  following  simple  estimate 
(26)  «uv-Uli  n  +  (rrjv)  1/2  IV*(uv-O),0fO  i  CI1-2V) 
where  I  Pi Q  refers  to  the  norm  in  the  quotient  space. 


0,n  ' 


Let  l  =  {T\}”_^  136  a  triangulation  of  01  with  the  property  that 
no  vertex  of  a  triangle  lies  in  the  interior  of  a  side  of  another. 

P (p) #0  denotes  the  set  of  continuous  functions,  whose  restriction 
to  each  triangle  is  a  polynomial  of  degree  at  most  p  and  whose  trace 


on  vanishes.  The  finite  dimensional  problems  we 


consider 


are  those  that  arise  when  the  energy  expression  (24)  is  minimized  in 
p(p],0  *p(p]»0  #  p  o.  For  values  of  v  away  from  H  the  finite 
dimensional  solutions  are  easily  seen  to  be  as  good  as  the  best  H^- 
approximations  to  u  from  within  the  spaces  xplPl'0  # 


The  aim  of  this  section  is  to  prove  an  optimal  asymptotic  error 
estimate  in  terms  of  negative  powers  of  p,  which  is  valid  in  the  limit  as 
v  approaches  %.  First  a  simple  approximation  result. 


Lemma  3.1.  For  any  k  >1  and  any  c  >  0  there  exists  a  constant 

C,  such  that 
k,  e 


inf  lv-Qpl.  <  C.  (p+l)“k+1+Elvt  _ 

qP6  p[P]'0*MP]'0  "  "  '  '  ' 


for  any  p  ^  0  and  any  v  €  (Hk(fl).  <">  • 


Proof :  see  [4]  . 

Let  u  denote  the  element  of  p^'®  that  minimizes 

-V,p  r  r 

(24) ,  we  shall  now  prove  that  u  converges  to  u  at  the  optimal 

— V  i P  —v 

rate  even  in  the  limit  as  v  approaches  >j.  Unlike  the  result  proven 
in  Ill),  for  the  case  of  a  smooth  boundary  an  and  natural  boundary 
conditions  on  all  of  an,  the  following  theorem  does  not  establish 
uniformly  valid  optimal  convergence  rates.  The  weaker  formulation 
is  due  to  the  nonsmoothness  of  the  boundary  an  .  for  a  fixe*' 

load  F  the  corner  singularities  and  therefore  the  smoothness  of  the 
solution  u^  is  likely  to  depend  on  v  .  This  means  that  there  are  no 
uniformly  valid  optimal  rates,  rather  the  optimal  rates  depend  on  v 
It  is  undoubtedly  possible  to  establish  error  estimates  that  exhibit 
these  optimal  rates  (depending  on  v  )  but  with  constants  that 
are  independent  of  v  .Technically  however  this  introduces  a  lot  of 
difficulties  that  are  irrelevant  to  the  main  scope  of  this  section, 
namely  to  understand  exactly  what  makes  the  p-version  superior  for 
nearly  incompressible  materials. 
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Theorem  3.1.  Assume  that  s  >  1  and  v .€  10, 1/2 [ -  Let  uv  be 
the  solution  to  (22)  and  (23) ,  and  let  U,  P  be  the  solution 
to  (25).  For  any  p  ^  0,  uy  denotes  the  element  in 
Pfp3,°  x  pi?1'0  that  minimizes  (24). 

Given  e  >  0  there  exists  C_  .  (independent  of  v  and  p)  such 
that  whenever  U  e  HS(fl)  » 

1 -v  ~  -v/Pl  i,o  +  (i^2v)1//2  ,v*<Hv"Sv,p)lo/n 
iCS/e<P+1»‘S+1+e  ‘“‘s^  +  C(1-2V)V2  ip1q ^ 

(the  constant  C  is  independent  of  all  parameters) . 


proof-  Let  T  (v.)  denote  the  projection  of  v  onto  P3p3  '®xp^p3 

-  v,p  —  — 

in  the  energy  inner  nroduct  associated  with  (24) .Then  u  _  =  T  (u  ) 

— v,P  v,P  -v 

and  therefore 

i  'Sv-U'i.n  * 

*  ,S‘tv,p,2>,i,b  *  *1^’  ,,',2'Tv.p®)'o.n 

*  ,Tv.P  ‘2-Vi.p*  lI^7>',,,-rev.p|2-Hv)i'o,0  • 


The  operators  Ty  are  uniformly  bounded  in  the  norms 

1  'l.n^  (FTv)1}  ,V‘(  3,ofn  and  by  using  the  estimate  (26) 
we  thus  conclude 
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•Sv-Sv.p'i.n  +  (l^v)1/2  l7*(^-iiv<p),o,n 
(27)  <  »g- Tv#p(U)*1>Q  +  (j^v)172  l*-<2  -Tv,p<S»>'o,fl 

+  C(l-2v)1/2  I  PI 0#n. 

Similar  to  Lemma  2.2  in  till  (for  which  we  did  not  use  the  smoothness 
of  the  domain  n)  we  get  for  any  v  with  V*v  =  0  that 

■I  -  ^.pCD'i.o  +  ‘irjv)1/2  ,7,(^  "  Tvfp(-,)'o,n 


<  D(Cp)  { 


qP6  plpl  »0  X  pfpl  »° 


•3-  sP'i.q)  ' 


2  1 

where  Cp  denotes  the  i(L  ,H  )-norm  of  any  maximal  right-inverse  for 
the  divergence  operator  acting  on  P ^  ' 0 x  P ^ ' 0  (here  we 
assume  thatp>3  so  that  V*  (Ptpl  ,0  *  Ptpl  ,0)  /  0)  .  Theorem  2.1 
together  with  Remark  2.  3  shows  that  we  can  obtain  C^<  C(p+l)K,and 
as  a  consequence  of  this  and  Lemma  3.1, 


<  CkjE(p+l)2K’k+1+E/2  |V|^  n  , 


k  °i 

provided  v  e  H  (£2)  n  H  (fl)  and  V*v  =  0  . 

The  operators  Tv  p  are  uniformly  bounded,  i.e.. 


(29)  p<3'i,0  *  irar'l’l'-<i-Tv,p(v))>0_I! 


i  c 


provided  V*v  ■  0. 
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Let  Hr  (£))  denote  the  space  Hr(G)n  {vjv-v  -  0)  ,  r  >  1  . 

Since  fl  is  simply  connected  Hr(ft)  may  also  be  characterized  as 
**fy  $»  -  «eHr  +  1(fi)  <"*  H2  tn>  i.e,  «*(())  is  isomorphic  to 

Hr+1(fi)n  h 2(fl).  Using  the  same  ideas  as  in  1 2]  (see  also  (12))  only 
with  the  biharmonic  operator  instead  of  the  Laplace  operator  one  gets 
the  following  result  concerning  complex  interpolation  and  the  spaces 
Hr+1(fi)  nfi2(fl) : 

(Hr+1(0)  nfi2(n),  H  t+1(fi)  n*H2(n))0  =  Hr+1+(  t-r)0(n)n  H2(fi) 

for  any  1  <  r  <  t  and  0  <  8  <  l  (even  though  fl  is  not  smooth) . 

Due  to  this  result  and  the  aforementioned  isomorphism  we  can  now 

interpolate  between  (28)  and  (29)  to  obtain 


'J-’v.pt)'!/ 


■)^l  V*(v-t 


—  v 


(v))l„  o 
rP  o,n 


(30) 


<  ck#^e(P+i)2K6-(k’1,e  + 


I  Vl 

-  l+(k-l)6,n 


for  any  0  <  6  <  1  and  any  v e H1+ (k_1) 6 (ft)  .  Choose  k  >  4K-^i  +  1 
and  define  0  =  .  The  estimate  (30)  with  U,  the  solution  to  (25), 

inserted  for  v  then  reads 

'SMWSMx,,,  *  liars’11  ,’-<E-*v,p«Sl»  1  o,ft 

(31) 


□ 


Combining  (27)  and  (31)  we  get  the  desired  estimate. 
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Remark  3.1.  An  alternate  way  to  formulate  the  result  of  Theorem  3.1 
is,  in  light  of  the  estimate  (26), 

_  — v,P* i,n  +  (rrs)  1/2  ,7-<E-!iv,p>VG 


<  c  e(p+l)‘s+1+e  IU» 


H's ,0  +  C(1-2V)1/2  IPl0ffl 


By  choosing  the  penalty  parameter  ^rjv  sufficiently  large  and 
minimizing  (24)  in  x  pfp),0  we  thus  obtain  a  solution  u„ 

-v,p 

that  converges  to  U,  the  solution  of  the  Stokes  problem,  at  the  rate 

»g4-l  g 

(p+1)  (modulo  an  e)  provided  U  €  h  (ft)  .  This  rate  is  optimal 


for  a  general  U  e  H  (ft)  ?  if  namely 

•a-Xj'x.o  ic,Pj+1)"8+l'4 

x  tPi],0  [p.),o 

for  some  <5  >  o  and  some  sequences  v ^  e  P  J  x  p  J  t 

Pj+1/Pj  5.  A  <  “•  Pj  +  "•  then  one  could  conclude  (cf.  [4])  that 

a 6  Ov 

for  any  t  <  6  and  any  triangle  e  Z.  (u  G  means  that 

U  €  Hs+t(K)  on  all  compact  subsets  K  of  7\.) 


Remark  3.2.  Whether  one  thinks  of  uv^  as  an  approximation  to  the 
displacement  for  a  nearly  incompressible  material  or  as  an  approxima¬ 
tion  to  the  solution  of  the  Stokes  problem,  the  behavior  of  the  error 
is  remarkably  different  from  that  of  more  standard  finite  element 
discretizations.  If  denotes  a  triangulation  of  mesh  size  h  and 
one  minimizes  (24)  in  the  set  of  continuous  piecewise  polynomials 
of  some  fixed  degree  for  decreasing  values  of  h,  then  it  is  well 
known  that  the  convergence  rates  obtained  this  way  are  not  optimal  as 
v  approaches  j,  for  piecewise  linear  functions  there  is  no  conver¬ 
gence  at  all.  We  refer  to  the  introduction  of  [11]  for  more  details. 


30 


References 

[1]  Babuska,  I.  and  Dorr,  M. ,  Error  estimates  for  the  combined  h  and 
p  versions  of  the  finite  element  method.  Num.  Math.  £7  (1981) 
pp.  257-277. 

[2]  Babuska,  I.,  Kellogg,  R.B.  and  Pitkaranta,  J.,  Direct  and  inverse 
error  estimates  for  finite  elements  with  mesh  refinements. 

Num.  Math.  33  (1979)  pp.  447-471. 

(3J  Babuska,  I.  and  Szabo,  B.A.,  On  the  rates  of  convergence  of  the 

finite  element  method.  Report  WU/CCM-80/2,  Washington  University, 
St.  Louis.  To  appear.  Intern.  J.  Num.  Meth.  Engng. 

[4]  Babuska,  I.,  Szabo,  B.A.  and  Katz,  I.N.,  The  p-version  of  the 

finite  element  method.  SIAM  0.  Num.  Anal.  18  (1981)  pp.  515-545. 

15]  Bergh,  J.  and  Lbfstrom,  J.,  Interpolation  Spaces.  Springer,  1976. 

[6]  Fix,  G.  J.,  Gunzburger,  M.D.  and  Nicolaides,  R.A.,  On  mixed 
finite  element  methods  for  first  order  elliptic  systems. 

Num.  Math.  3[7  (1981)  pp.  29-48. 

[7]  Katz,  I.N.,  Peano,  A.G.  and  Rossow,  M.P.,  Nodal  variables  for 
complete  conforming  finite  elements  of  arbitrary  polynomial  order. 
Comp.  6  Maths,  with  Appls.  £  (1978)  pp.  85-112. 

18]  Morgan,  J.  and  Scott,  R.,  A  nodal  basis  for  piecewise 

polynomials  of  degree  >  5.  Math.  Comp.  29  (1975)  pp.  736-740. 

[9]  Raviart,  P.A.  and  Thomas,  J.M.,  A  mixed  finite  element  method 
for  2-nd  order  elliptic  problems,  pp.  292-315  in  "Mathematical 
Aspects  of  Finite  Element  Methods".  (Galligani,  I.  and  Magenes,  E. , 
Eds.)  Springer  Lecture  Notes  in  Math.  606,  1975. 

110]  Temam,  R. ,  Navier-Stokes  Equations.  North-Holland,  1977. 


3 


111]  Vogelius,  M. ,  An  analysis  of  the  p-version  of  the  finite 

element  method  for  nearly  incompressible  materials.  Uniformly 
valid,  optimal  error  estimates.  (To  appear,  Num.  Math.) 

(12)  Vogelius,  M.  and  -  Papanicolaou  ,  fi. ,  A  projection  method  applied 
to  diffusion  in  a  periodic  structure.  (to  appear,  SIAM  J. 
Appl.  Math.) 


SECURITY  CLASSIFICATION  OF  THIS  RAGE  (Whon  Data  BntaraB) 


REPORT  DOCUMENTATION  PAGE 


4.  TITLE  (and  Subtitle) 

A  RIGHT- INVERSE  FOR  THE  DIVERGENCE  OPERATOR  IN 
SPACES  OF  PIECEWISE  POLYNOMIALS.  APPLICATION 
TO  THE  p- VERSION  OF  THE  FINITE  ELEMENT  METHOD. 


7.  AUTHORf*; 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


1.  RECIPIENT'S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  A  PERIOD  COVERED 

Summary  Report  -  no  specific 
reporting  period 


<>  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERfaj 


Michael  Vogelius 


r  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Mathematics  Research  Center,  University  of 
610  Walnut  Street  Wisconsin 

Madison.  Wisconsin  53706 


5  I  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


(See  Item  18  below) 


DAAG29-80-C-0041 

N00014-81-K0002 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Work  Unit  Number  3  - 
Numerical  Analysis  and 
Computer  Science 


12.  REPORT  DATE 

May  1982 


IS-  NUMBER  OF  PAGES 


MONITORING  AGENCY  NAME  A  ADORESSf"  tUHoront  bom  Controlling  Ollleo)  IS.  SECURITY  CLASS,  (of  thlo  report; 

UNCLASSIFIED 


IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


IS.  DISTRIBUTION  STATEMENT  (o I  thlo  Report) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (o I  tho  eft.fr eel  entered  In  Block  30,  II  different  from  Roport) 


IS.  SUPPLEMENTARY  NOTES 

U.  S.  Army  Research  Office  < 

P.  O.  Box  12211  i 

Research  Triangle  Park 
North  Carolina  27709 


18.  KEY  WORDS  ( Conttnu •  on  rtv*rM  •/<#•  It  no coaaaay  and  Identify  by  block  nuatbor) 


Office  of  Naval  Research 
Arlington,  VA  22217 


discretization  of  PDE's,  finite  element  method,  piecewise  polynomial  spaces, 
divergence  operators,  optimal  convergence  rates,  nearly  incompressible 
materials 


20.  ABSTRACT  (Continue  on  roooroo  •!«  If  nacmaaary  and  Identity  by  bloc*  number) 

'"’‘in  the  first  part  of  this  paper  we  study  in  detail  the  properties  of 
the  divergence  operator  acting  on  continuous  piecewise  polynomials;  more 
specifically,  we  characterize  the  range  and  prove  the  existence  of  a  maximal 
right-inverse  whose  norm  grows  at  most  algebraically  with  the  degree  of  the 
piecewise  polynomials.  _  _ - 


1473  COITION  OF  I  NOV  SS  It  OBSOLETE 


^  (cont.) 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWRaa  Data  Bntorod) 


fcfrr .  .  '|jji 


ABSTRACT  (cont.) 

n  the  last  part  of  this  paper  we  apply  these  results  to  the  p-version 
of  the  Finite  Element  Method  for  a  nearly  incompressible  material  with 
homogeneous  Dirichlet  boundary  conditions.  We  show  that  the  p-version 
maintains  optimal  convergence  rates  in  the  limit  as  the  Poisson  ratio 
approaches  1/2.  This  fact  eliminates  the  need  for  any  'reduced  integration* 
such  as  customarily  used  in  connection  with  the  more  standard  h-version  of 
the  Finite  Element  Method 


